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Abstract. We report on a study of the 7-ray continuum 
emission from the Crab supernova nebula and on a search 
for nuclear de-excitation 7-ray lines. Crab is the brightest 
continuum source in the 1-10 MeV 7-ray sky, and its con- 
tinuum radiation is most likely of synchrotron origin. It is 
a likely source of cosmic rays through shock acceleration 
and thus a potential candidate for 7-ray line emission from 
nuclear interactions. Five years of COMPTEL observations 
enable a fine spectral binning to investigate the behaviour of 
the 0.75-30 MeV emission in detail and to search for nuclear 
de-excitation lines on top of the continuum. The nebular 
spectrum shows a break at the edge of the COMPTEL energy 
range and connects well to the EGRET spectrum, probably 
reflecting electron energy losses in the synchrotron emission 
scenario. Such a smooth continuum model alone may not be 
sufficient to explain the observations. A weak bump in the 
spectrum at 1-2 MeV may be present. No significant evidence 
for distinct line emission is seen, but the presence of a blend 
of line features or another synchrotron component cannot be 
excluded. 
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1. Introduction 

Strong shocks in Galactic supernova remnants (SNRs) are 
widely believed to be sources of cosmic rays (CRs) below 10*^ 
eV. This has only been substantiated, however, through de- 
tection of accelerated energetic electrons by their synchrotron 
emission. Direct evidence for the acceleration of the nuclear 
component is still to be found. Searches for signatures of the 
nuclear component in the TeV range (e.g. Atoyan & Aharonian 



also be deduced from the presence of 7-ray lines which these 
particles produce when interacting with the ambient medium. 

Nuclear interaction line emission has so far been reported 
from solar flares (e.g. Chupp 1990; Murphy et al. 1991 ; Ryan et 
al. 1992) and in only one occasion from an extra-solar source. 



the Orion Complex (Bloemen et al. 1994a, 1997c). Evidence is 
seen now in emission from the inner Galaxy as well (Bloemen 
et al. 1997b)| ). These detections were obtained with the Comp- 
ton Telescope (COMPTEL) aboard the Compton Gamma Ray 
Observatory (CGRO). COMPTEL operates between 0.75 and 
30 MeV, which is the energy regime where signs of nuclear in- 
teractions can be expected. The observed 3-7 MeV emission 
from Orion was tentatively attributed to ^^C* (4.44 MeV) and 
*^0* (6.13 MeV) de-excitation lines. Solar flare spectra show 
several additional lines from heavier nuclei (^''Ne, ^*Mg, ^^Si), 
mainly in the 1-2 MeV regime. Here we report on a search for 
such 7-ray lines from the Crab nebula, using COMPTEL data 
obtained during the first five years of the mission. 

The Crab nebula is an extensively studied 943 year old 
SNR. The pulsar wind from the central, isolated pulsar sup- 
plies high-energy particles that synchrotron radiate (Shklovsky 



1953; Dombrovsky 1954; Oort & Walraven 1956) while gy 



rating in the magnetic field of the nebula (Kennel & Coroniti 
1984). However, in situ acceleration has to occur as well. The 



1996; Allen et al. 1995), for instance, have only provided up- 



per limits. The existence of accelerated protons and nuclei can 



radio-to-7-ray spectrum consists of a mix of differe nt pow er- 
law components, separated by breaks (e.g. Zombeck 1982 ; De 
Jager et al. 1996| ), which seem to indicate three electron popu- 
lations. The continuum emission up to ~ 1.8 x lO^'^Hz is as- 
sociated with the optically visible nebula. The break may re- 
flect the energy losses for electrons as old as the SNR; the av- 
erage field of the optical nebula can then be calculated to be 
- 3 X 10~4 G (e.g. Trimble |1982| ; Marsden et al. |1984| ). Conse- 
quently, the electrons responsible for the torus of X-ray emis- 
sion, the spectrum of wh ich breaks somewhere between 60 and 
150 keV (Bartlett |1994 . should have a lifetime of ~ 1 month. 
The inner edge of the nebula lies at 8" (0.08 pc), so it is clear 
that these electrons have been accelerated in or near the nebula. 
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De Jager et al. (1996) discuss evidence for a third spectral 
break at ^ 26 MeV, from early observations by COMPTEL and 



2 



R.D. van der Meulen et al.: COMPTEL 7-ray study of the Crab nebula 



Table 1. COMPTEL observations used in the present Crab 
analysis. The observation numbers are given in standard CGRO 
notation. 



Obs 




Tend 


Obs 




Tend 


# 


dd-mm-yy 


dd-mm-yy 


# 


dd-mm-yy 


dd-mm-yy 


ga 


28-04-91 


07-05-91 


321.5 


15-02-94 


17-02-94 


1.0 


16-05-91 


30-05-91 


337 


09-08-94 


29-08-94 


31 


11-06-92 


25-06-92 


412 


28-02-95 


07-03-95 


36.0 


11-08-92 


12-08-92 


413 


07-03-95 


21-03-95 


36.5 


12-08-92 


20-08-92 


419.1 


04-04-95 


1 1-04-95 


39.0 


01-09-92 


17-09-92 


419.5 


09-05-95 


23-05-95 


213 


23-03-93 


29-03-93 


420 


23-05-95 


06-06-95 


221 


13-05-93 


24-05-93 


426 


08-08-95 


22-08-95 


310 


01-12-93 


13-12-93 


502 


17-10-95 


31-10-95 


321.1 


08-02-94 


15-02-94 








" Data obtained during 


the verification phase 







EGRET (also aboard CGRO), which is addressed below. This 
would indicate another population of energetic electrons, with 
lifetimes of 1 day. We present further evidence for this spec- 
tral break using more COMPTEL and EGRET observations. 

2. Instrumentation, observations & data analysis 

2.1. Instrument 

COMPTEL is designed to detect 7-ray photons with energies 
between about 0.75 and 30 MeV with an energy resolution of 
5-10% FWHM. It has a field-of-view of 1 steradian and 
an angular resolution of typically l°-3°, which enables the 
instrument to monitor several objects simultaneously. The lo- 
cation accuracy is about 0.5° for a strong source. A detailed 
description of the detection principle and instrument is given 
by Schonfelder et al. ( 1993 ). Ideally, incoming 7-ray photons 
are first Compton scattered in an upper detector layer and then 
completely absorbed in a lower detector layer. The measured 
energy deposits and locations in these layers determine the 
scatter direction, scatter angle, and total energy of each photon. 
For selected energy bands, the telescope events are binned in a 
3-dimensional data space, which consists of two scatter direc- 
tion coordinates and a Compton scatter angle coordinate. Here 
we use a binning ofl° x 1° x 1°. The source response function 
of the instrument in this data cube has a cone-like shape and 
depends on the actual source spectrum. In the present work we 
have adopted an E^"^ power-law input spectrum, but our find- 
ings are not sensitive to this specific choice. 

We applied a maximum-likelihood method (de Boer et al. 



1992) to obtain flux estimates and to construct spectra with nar- 
row energy bins equal to twice the FWHM of the COMPTEL 
energy resolution, given by 



FWHM(£:) = 0.0236 x (14.61£; + 2.biE^) 



•2^,1/2 



with E in MeV (Schonfelder et al. 1993). Our method of an- 
alysis does not allow a significantly finer energy binning at 
this stage. At high energies (> 10 MeV) the limited number 
of events requires even larger bins. In this study we have a total 
of 30 statistically independent bins. 



Table 2. Crab nebula (as 


plotted in Fij 


5. 1) and total (Fig. 


fluxes. 










E 


FluXunn 


Errorunp 


Fluxtot 


Errortot 


MeV 




10-^7 


— 2 —1 

cm s 


0.78-0.96 


34.07 


1.57 


42.75 


1.00 


0.96-1.16 


31.11 


1.30 


36.68 


0.82 


1.16-1.38 


21.54 


1.09 


27.80 


0.69 


1.38-1.62 


18.68 


1.04 


23.78 


0.65 


1.62-1.88 


13.03 


0.82 


17.02 


0.52 


1.88-2.16 


8.90 


0.71 


11.62 


0.45 


2.16-2.48 


7.84 


0.70 


10.25 


0.44 


2.48-2.84 


8.66 


0.57 


9.79 


0.36 


2.84-3.22 


6.14 


0.52 


7.28 


0.33 


3.22-3.62 


4.41 


0.46 


6.13 


0.29 


3.62-4.08 


4.28 


0.41 


5.80 


0.26 


4.08^.56 


4.44 


0.36 


4.77 


0.22 


4.56-5.08 


2.94 


0.33 


3.96 


0.21 


5.08-5.66 


3.42 


0.32 


3.69 


0.20 


5.66-6.26 


2.30 


0.28 


2.91 


0.18 


6.26-6.94 


2.15 


0.28 


2.64 


0.18 


6.94-7.64 


1.60 


0.25 


2.11 


0.16 


7.64-8.42 


1.17 


0.20 


1.80 


0.13 


8.42-9.26 


1.36 


0.17 


1.75 


0.11 


9.26-10.16 


1.24 


0.15 


1.47 


0.10 


10.00-11.20 


1.05 


0.15 


1.57 


0.10 


11.20-12.48 


1.12 


0.14 


1.37 


0.09 


12.48-13.92 


1.12 


0.14 


1.24 


0.09 


13.92-15.52 


0.62 


0.13 


0.92 


0.09 


15.52-17.28 


0.83 


0.14 


0.96 


0.09 


17.28-19.28 


0.50 


0.14 


0.68 


0.09 


19.28-21.60 


0.31 


0.14 


0.69 


0.10 


21.60-24.08 


0.23 


0.15 


0.29 


0.10 


24.08-26.88 


0.16 


0.18 


0.26 


0.11 


26.88-30.00 


0.13 


0.27 


0.09 


0.18 



On average, more than 95% of the number of events tagged 
by the instrument as Compton-scattered 7-ray should be at- 
tributed to the background (instrumental and isotropic). The 
likelihood analysis requires, for each energy interval sepa- 
rately, a careful estimate of this background for which a variety 
of methods has been studied. The results shown in this paper 
are based on a method of background determination in which 
a filter technique is applied to the data space. This method is 



described by Bloemen et al. (1994b), although we applied here 



an improved algorithm involving an iterative process of back- 
ground estimation. The filter eliminates to first order any source 
signature present. In each iteration the background is further 
corrected for the smeared-out source signature. 

In this work we use point spread functions (PSFs) from an- 
alytical modeling based on single-detector calibrations. PSFs 
from Monte Carlo simulations of the instrument are preferred, 
but not available yet for the narrow energy bins used here. A 
globally somewhat softer spectrum is expected with the simu- 
lated response (about 0.1 in spectral index). 
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Fig. 1. COMPTEL s pectru m of the Crab nebula, together with GRIS low-energy data (Bartlett et al. |1993| ) and EGRET high- 
energy data (Fierro |1996| ). Up to 10 MeV, the bin width of the COMPTEL data points is equal to twice the FWHM of the 
instrumental energy resolution. Between 10 and 30 MeV, broader energy bins were chosen (Sect. 2.1). The two COMPTEL 
upper limits at the highest energies (Table 2) are not shown. 



2.2. Observations 



3. Results 



Crab was within 30°of the instrument pointing during several 
observations in the first 5 years of CGRO operations (Table |l]). 
Observation 2.5 was excluded because the instrumental settings 
were optimized for solar flare observations. 



2.3. Phase selection 

In addition to the Crab nebula, the Crab pulsar is a strong 
(pulsed) 7-ray source as well, which fully dominates the high- 



Nolan et al. 

range, about 20% of the total Crab emission is estimated to 



energy (> 100 MeV) 7-ray emission (e.g. Kanbach et al. 1990; 

[199^ Fierro [l996| ). In the COMPTEL energy 



be due to the pulsar (cf Fig. 6 in Much et al. 1995a). COMP 



TEL cannot resolve the Crab pulsar from the surrounding neb- 
ula spatially, but it is possible to disentangle them by pulsar 
phase selection, assuming that there is no pulsar emission in 
the "off-phase at 7-ray energies. We determined the unpulsed 
emission by selecting the events in th e off-pulse phase 0.525- 
0.915, as defined by Nolan et al. ( I1993I ). Our analysis thus gives 
39% of the flux, which is then normalized to the full period. 
The observations listed in Table |l] were combined and folded 
by our pulsar analysis software. 



Figure 1 shows a power-per-decade spectrum of the off-pulse 
fluxes (Tab. 2). Our data points are in good agreement with the 
COMPTEL spectrum presented by Much et al. ( |l995b| ), which 
was based on observations obtained during the first two years 
of the mission and contains 7 data points only. Statistical error 
bars are shown. The absolute calibration uncertainty is conser- 
vatively estimated to be 30%. The total number of photons in 
this s pectru m is ^ 68000. Measurements from GRIS (Bartlett 



1993) and EGRET (Fierro 1996, April 1991 - Aug. 1994) 



et al. 

are added to place the COMPTEL points in a broader perspec- 
tive. 



De Jager et al. ( 1996 ) have modelled the emission of the 
Crab nebula. They used an inverse-Compton (IC) component 
and a synchrotron comp onent w ith an exponential cut-off to fit 
COMPTEL (Much et al. |l995b| ) and EGRET data (both April 
1991 -May 1993): 

dN/dE = Ks{E/3.5MeY)-'^^exp{-E/Eo) 



-/s:ic(£;/1000MeV) 



-Fic 



Their best-fit result is shown as a dashed line in Fig. 1 (lim- 
ited to the 1-100 MeV range). Figure 1 shows that the spectral 
break represented in the model begins in the COMPTEL do- 
main. We have fit our data and the first 3 EGRET points to the 
first term of the equation. We ignore the IC component of the 
model, which is valid because the IC influence on the used data 



4 



R.D. van der Meulen et al.: COMPTEL 7-ray study of the Crab nebula 



.001 

.003 



.001 

.003 



> 

(U 



o 
1) 

X 
3 



.001 
.003 



.001 
.003 



"■"■"1 I 

a) Crab off-pulse 



H — h 



H 1 1 1 ' t I 



b) Crab off-pulse 



H — h 



+ , + + 

— f I "I" I r 



H h 



c) Simulation Crab off— pulse 



+ + 



H — h 



H 1 — I— I- 



d) Simulation Crab off— pulse 



H — 1— 1- 



H — h 



H 1 — I— I- 



e) Simulation Crab off— pulse 



Energy (MeV) 

Fig. 2. a) 0.75-9.7 MeV COMPTmi spectrum of the Crab neb- 
ula (solid lines). The bin sizes are twice the FWHM of the 
instrumental energy resolution and shifted by 1 FWHM with 
respect to the points in (b), also shown (dotted), b) 0.78-10.2 
MeV nebula spectrum as in Fig. 1; the triangles indicate posi- 
tions of some candidate nuclear interaction lines : ^"'Mg (1.37 
MeV), 20Ne (1.63 MeV), ^^Si (1. 78 Me V), ^^C (4.44 MeV), 
and (6.13 MeV) (Ramaty et al. |1979|) c,d,e) Simulations of 



E ^ spectra at a level of 1.38 X 10 MeV cm ^ s ^. 



points is negligible. For the fit parameters we find: the energy 
cut-off £;o = 41 ± 3 MeV, Ks = (1.29 ± 0.04) x lO"'' cm-^ 
s^i MeV"\ and = 2.02 ± 0.03, with a reduced ranging 
from 1.9 to 2.2 for 28 degrees of freedom (dof). The latter 
value was obtained by fitting the model to a spectrum with bins 
shifted by 1 FWHM of the energy resolution. Figure 2a shows, 
in solid lines, this shifted spectrum from 0.75-10 MeV. For 
comparison, the unshifted spectrum, also shown in Fig. 2b, is 
represented by dotted lines. The Icr errors on the fit parameters 
have been calculated by increasing the Xmin ^i'^^ 3.5 (Lamp- 
ton et al. 1976 ). The result can be seen in Fig. 1 (solid line). 
The cut-off energy is somewhat sensitive to the EGRET data 



points. Leaving all EGRET points out gives fit parameters of 
Ea = 20t^ MeV, A'^ = (1-47 ± 0.04) x lO"'^ cm^^ g-i 
MeV"\ and L^ = 1.90 ± 0.03, with a reduced fit rang- 
ing from 2.0 to 2.3, for 25 dof. For comparison, De Jager et al. 



(1996), who included the COMPTEL calibration uncertainty 



in their fit but excluded the < 1 MeV COMPTEL point, found: 



MeV-\ and 



Eq = 26tf MeV, Ks = 1.25 x 10" - cm 
Ls = 1.74 ±0.42. 

The fit is not perfect, as indicated by the x^ value obtained 
above, but no obvious systematic trends can be seen. The devi- 
ations from the model fit are below the 3 a level in the individ- 
ual energy bins. In order to assess the systematic errors of our 
method we analyzed three Monte Carlo simulations of the Crab 
up to 10 MeV (Fig. 2cde), assuming an E^^ spectrum with an 
integrated 1-10 MeV flux of 1.38 x 10~^ cm^^ s^^, as derived 
from the observations. It appears that with this input flux (cor- 
rected for the phase selection, i.e. 0.39 x 1.38 x lO^'' cm^^ 
s^^) our method generates randomly placed features similar to 
those in the observed spectrum. Therefore, no significance can 
be given to the observed features. Best fits of the simulations to 
power-law spectra 

dN/dE = K{E/3.5MeYy'^ 

give Xrod values of 0.9, 1.4, and 0.9 respectively, for 18 dof. 
For comparison, the best power-law fit to the observed data up 
to 10 MeV has an index of 2.02^^;!^^, K ==(1.12 ±0.03) x 10"'' 
cm^^ s^^, and a x^d ranging from 2.6 to 3.0, depending on the 
bin positions. The la errors on the fit parameters have been cal- 
culated by increasing the Xmin ^i'^^ 2.3 (Lampton et al. 1976). 



Since a mixture of line features (including broad-line com- 
ponents) may well be present, it cannot be excluded that we 
overestimated the continuum contribution. On the other hand, 
a decrease of the fitted continuum level would not exhibit clear 
evidence for known candidate lines (triangles in Fig. 2b) with- 
out introducing other features as well (e.g. near 2.8 MeV). 
Gravitationally redshifted 7-ray lines caused by ions falling 
onto the pulsar have not been included in the candidate lines, 
because it is unknown how deep in the atmosphere the particles 
would interact, a nd thu s how large the redshift would be. 

Much et al. (1995a) reported that the Crab unpulsed emis- 
sion might be time variable. Spectral changes (e.g. the position 
of the break energy) may occur on short time scales, possibly 
affecting the global appearance of our spectrum. Time variabil- 
ity of the Crab emission will be the topic of future study. 

In order to obtain the cleanest data, we selected only 39% 
(the unpulsed fraction) of the available nebula events. How- 
ever, as the emission from the nebula dominates the total emis- 
sion from the Crab, we have used the set of all available (pul- 
sar contaminated) events for our analysis as well. Figure 3d 
shows the total Crab spectrum using all observations; Fig. 3a 
shows the same spectrum (dotted lines) and the spectrum with 
shifted bins (solid lines). The integrated 1-10 MeV flux of 
1.75 X 10^^ cm^^ s^^ has been used as input parameter for 
simulations of an E^^ spectrum (Fig. 3efg). The best power- 
law fits to the simulations have Xrod values of 1.5, 2.4 and 
1.7 (18 dof). Clearly, this is not perfect; the deviations may 
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result from the fact that we have used the actual observations 
as the basis for our simulations, so that the smoothed Crab sig- 
nal may influence our findings. Anyway, the best power-law 
fit to flie observed data up to 10 MeV [T = 2.02 ± 0.01, 
K =(1.41 ± 0.02)xl0"'^ cm"2 s^^], has a x^cd ranging from 
5.0 to 7.4 (depending on the bin positions), which is a clear in- 
dication that the spectrum is more complicated. The main dif- 
ference between observation and simulations is in the low en- 
ergy range; at 1-2 MeV there is a hint for a broad feature. Some 
evidence for this feature could in fact already be seen in the un- 
pulsed spectrum of Fig. 1 , but it seems more pronounced in the 
total Crab spectrum. After splitting the data in two time inter- 
vals (Fig. 3bc, observations 0-337 and 412-502 respectively), 
some evidence for this feature is seen in both spectra. We can- 
not exclude, however, that it is associated with the pulsar, if 
real. The B ATSE (35-1700 keV) Cycle 1-3 total Crab spec- 
trum (Ling |1997| ), altho ugh h igher than the ORIS and COMP- 
TEL points (Much et al. 1996 ), shows a deviation from a simple 
power-law spectrum, which can be interpreted as a rise near 1 
MeV as well. If real, the feature could well be due to a blend 
of lines, or it could be the synchrotron signature of a different 
population of electrons. Morphological comparisons of OSSE, 
BATSE and COMPTEL Crab spectra will be the topic of future 
study. 

4. Conclusions 

COMPTEL observations from the first five years of the mission 
were combined to search, in a fine binned spectrum, for nuclear 
de-excitation 7-ray lines from the Crab supernova nebula. No 
significant evidence for line emission is found. A reasonable 
fit to the off -pulse fluxes can be obtained with the sy nchrotron 
component of the model from De Jager et al. (1996 1, giving a 
break energy of about 25^0 MeV. An indication for a bump 
at 1-2 MeV cannot readily be explained by known systematic 
and statistical errors. It is possible that line contributions or an 
extra continuum component have resulted in an overestimation 
of the single continuum component. 
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Fig. 3. a) 0.75-9.7 MeV COMPTEL spectrum of the total Crab 

emission (solid lines). The bin sizes are twice the FWHM of the 
instrumental energy resolution and shifted by 1 FWHM with 
respect to the points in (d), also shown (dotted). b,c,d) 0.78- 
10.2 MeV total Crab spectrum using observations from respec- 
tively the first three years, the fourth and fifth, and the first five 
years of operations. e,f,g) Simulations of spectra at a level 
of 1.75 X 10-3 MeV cm'^ s'^. 



